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This paper reviews a recent non-Landau- Girizburg- Wilson (LGW) approach to super- 
fluid to Mott insulator transitions in two dimensional bosonic lattice systems, using a dual 
vortex field theory. 1 '' 2 ' The physical interpretation of conventional LGW theory of quan- 
tum criticality is re-examined and similarities and differences with the vortex picture are 
discussed. The "unification" of various competing (insulating) orders, and the coincidence 
of these orders with the Mott transition are readily understood in this formulation. Some 
aspects of the recent theory of "deconfined" quantum criticality, which are to an extent sub- 
sumed in this approach, are discussed. A pedagogical presentation of the "nuts and bolts" 
of boson-vortex duality at the hamiltonian level is included, tailored to a condensed matter 
audience. 

§1. Introduction 

Condensed matter theory continues to be challenged by the remarkable behavior 
of strongly correlated materials. The archetype of strong correlation physics is the 
Mott insulator, in which insulating behavior arises not due to the presence of filled 
bands (ultimately the Pauli exclusion principle), but to charge localization by strong 
local Coulomb interactions. Locally this physics is extremely simple, but extending it 
away from the atomic context, to the many-body problem is not so straightforward. 
The theoretical treatment of Mott insulators - and systems proximate to a Mott state 
- is complicated by the fact that such Mott localization reflects neither a feature 
of some electron-like quasiparticle spectrum, nor any kind of symmetry breaking. 
Thus the two workhorses of solid state physics, Landau's Fermi liquid theory and 
the symmetry/order parameter description of phases of matter, are not helpful in 
describing the basic Mott physics. 

In practice, most Mott insulators order at low temperatures, either magneti- 
cally, or by charge-order (charge density wave, stripe, etc.). These various orderings 
can be viewed as "competing" with one another, and with electronically-extended 
phases in which the Mott localization physics is not operative. Such competing or- 
ders are increasingly observable at the microscale through improvements in STM 
microscopy, 3 -^ 7 ) and in crystal and surface quality. A complication for theoretical 
treatments of competing orders is that there are too many different charge and/or 
spin ordering patterns, often of considerable complexity, consistent with the sim- 
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pie local Mott physics requiring only single occupancy of some orbitals. Therefore 
we view these orders as derivative phenomena, occuring as a consequence of Mott 
localization, not vice-versa. 

Nevertheless, competing orders seem to be an essential accompaniment to Mott 
localization. Indeed, there have been a number of "proofs" of late (some of which 
may even be rigorous!) that seem to require order of some sort in a Mott state which 
does not have an average electron occupation which is even per unit cell. 8 )' 9 ) More 
specifically, these proofs require ground state degeneracies in large systems, which 
may be associated physically either to broken symmetry (presumably the usual case) 
or to more exotic "topological order" . Thus a minimal requirement of any reasonable 
theory of the Mott transition is that it should lead automatically to conventional or 
at least topological order. 

In this paper, we discuss the Super Fluid (SF) to Mott transition in two-dimensional 
lattice boson systems. 10 ) This problem has recently become a very active experimen- 
tal area with the advent of cold trapped atoms in an optical lattice. 11 ) To connect it 
theoretically to the electronic Mott physics discussed above, one can view the bosons 
as tightly-bound "Cooper pairs" ; the above arguments for the necessity of competing 
order in the insulator thus apply when the boson filling non-integral. We give a some- 
what pedagogical review of recent work 1 )' 2 ) which satisfies the above requirement of 
naturally describing the emergence of competing orders in the insulator, while not 
putting this in "by hand". We focus on rational fillings / = p/q, with p, q relatively 
prime, and q > 1. This is accomplished by the use of duality, 12 )~ 14 ) a technical 
transformation of the hamiltonian which has the utility of allowing one to approach 
the SF-Mott (putative) quantum critical point from the superfiuid side, the opposite 
of the Landau-Ginzburg- Wilson (LGW) approach to superfluidity, 10 ) which expands 
about the normal phase in terms of the superfiuid order parameter. We argue that, 
because of the order in the insulator, such an LGW approach is unnatural. 

Thinking of the Mott transition as a Quantum Critical Point (QCP), understand- 
ing it is just the search for an appropriate continuum quantum field theory. Quantum 
field theories are "second quantized" descriptions of particles: point-like excitations 
created and annihilated by the quantum fields. LGW theory takes, loosely speaking, 
these particles to be the original lattice bosons. The dual approach described in 
this paper uses instead the vortices of the superfiuid as these "particles" . Because 
it is formulated in terms of vortex excitations, this approach takes advantage of the 
fact that the nature of the superfiuid ground state (in particular its symmetries) is 
independent of the boson filling. 

What begins as a worry - the non-local nature of the phase gradient/superflow 
surrounding a vortex - ends up as an advantage in the approach. In the vortex 
formulation, this non- locality is completely accounted for by a dual gauge field A^, 
which couples minimally to the vortices. The gauge field accounts for the phase 
winding of the boson wavefunction upon encircling a vortex. Turned around, the 
vortex wavefunction must wind when encircling a boson. The phase winding due 
to the average boson number encircled by a minimal motion of a vortex on the 
lattice, is captured in the dual theory by an Aharonov-Bohm "flux" of 2irf in 
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through a dual plaquette.*) When the bosons are at non-integer filling, this modifies 
the low-energy vortex dynamics in an essential way. In particular, it forces the 
vortices to appear in multiplets of q flavors, described by a vector of vortex fields, 
ipi, £ = 0...q — 1. These multiplets transform under a projective representation 
of the physical lattice symmetry group, a generalization of the ordinary notion of a 
representation, in which the group multiplication table is obeyed only up to phase 
factors, i.e. if a phase is "projected out". An important example is the x and y 
lattice translations on the square lattice, which obey 

T x T y = LoT y T x , (1-1) 

with to = e 2m f . This mathematically captures the Aharonov-Bohm phase described 
above. The modified group is dubbed a Projective Symmetry Group, or PSG. 18 ) A 
PSG is possible only because of the gauge nature of the vortex theory, "projection" 
being allowed by the lack of independent physical meaning of the local phase of the 
vortex field. 

The PSG dictates the form of the dual effective action, similarly to how the 
ordinary symmetry group dictates the free energy in LGW theory. One finds (see 
Sec© 

S = j d 2 rdT^{e^ x d v A x f + Y, [\{d» - iA»)<p t \ 2 + f\<p t \ 2 ] +£ in t}, (1-2) 

where C\ n t represents quartic and higher order terms in the {<pe}- , which are strongly 
constrained by the PSG. The form of C iri t is, however, specific to each value of q (see 
Ref. 1) for a general discussion, and Sec. ED for examples). 

Eq. (|1-2|) describes the superfluid/Mott physics through the gauge field A^. For 
instance, the phason mode of the superfluid is described in Eq. 1)1 -2jl as the gapless 
transverse "photon" mode on the gauge field. The condensation of any of the <pi 
fields leads to a "Higgs" mass for the gauge field, corresponding to the loss of the 
photon and the gap in the Mott phase. Competing "charge" orders are less apparent 
in Eq. (|l-2j) . but they are in fact encoded in the structure of the PSG. 

In particular, the order parameters for the different charge ordering patterns 
occuring in Mott states can be written in terms of "density wave" amplitudes, pq, 
describing the (complex) amplitude of a plane-wave oscillation in the charge density 
at wavevector Q. The non-trivial vortex PSG provides a link between the vortex 
multiplet and spatial symmetry operations. In fact, one can explicitly construct a 
set of such density wave operators, with 

Q mn = 2irf(m,n), (1-3) 

where m, n are integers. In particular, we find (see Sec. I5.3|) 

3-1 

Pmn = P Qmn = S (|Q mn |) u/™/ 2 £ <p* m+n J m . (14) 
e=o 

*' The same approach has been expored in Refs. 15)— 17). 
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In associating the p mn with a density, there is a general 'form- factor', S(Q), which 
cannot be determined from symmetry considerations, and has a smooth Q depen- 
dence determined by microscopic details and the precise definition of the density 
operator. 

Unlike in Landau theory, the density-wave order parameters describing the pos- 
sible "competing orders" in the Mott state are quadratic rather than linear in the 
quantum fields of the theory. The ordering is thus "weaker" than might be expected 
near a LGW-type charge ordering transition, consistent with the notion that the 
vortex action describes the Mott transition first, and competing charge orders in the 
insulator only as a secondary consequence. Describing both phenomena - the loss 
of superfluidity and onset of charge order - simultaneously is already an achieve- 
ment. It is possible because, since the ip# field carries the dual gauge charge, vortex 
condensation can describe both Mott physics (through the Higgs mass of when 
((fi) 7^ 0) and the emergent competing order arising through the p mn . An intriguing 
conseqence is that, even in the superfluid phase, one expects to see density wave 
modulations appearing in the vicinity of a localized vortex. ^ Equivalent but com- 
plementary viewpoints of this non-LGW quantum criticality for the special case of 
half-filling have been extensively developed in recent work on "deconfined quantum 
criticality" 19 ^ (see also Sec. 16.1.1)) . 

This paper is written to provide a "gentler" introduction to the work of Ref. 1). 
To keep it pedagogical, it was not possible to go beyond this work to discuss the 
intended applications of these theoretical ideas to electronic systems near a super- 
conductor to Mott insulator transition, with an eye to the under-doped cuprate 
materials. We cannot resist, however, pointing the reader toward this interesting 
direction. In Ref. 2), it was demonstrated that the same dual critical field the- 
ory applies to a somewhat more microscopically faithful representation of strong 
electronic pairing, the doped quantum dimer model. 20 ) The authors believe that it 
indeed applies more generally to any two-dimensional clean singlet superconductor to 
Mott insulator transition, in which the gap in the superconducting state is complete. 
Most recently, it was proposed that the observation of "checkerboard" charge corre- 
lations near vortices in BSCCO by Hoffman et ai, 21 ) is indicative of similar physics. 
This suggestion was used to extract some bounds on the inertial mass of a vortex 
from these experiments. 22 - 1 The reader should note that, because the cuprate ma- 
terials are gapless, d-wave superconductors, this interpretation goes boldly beyond 
the existing theory. Theoretical work to directly address the additional quasiparticle 
physics in gapless superconductors is ongoing. 

The remainder of this paper is organized as follows. In Sec. |2J we describe the 
general structure of simple models of lattice bosons, which contain superfluid to 
Mott transitions. Sec. describes the conventional LGW theory of the SF-Mott 
QCP for integer boson filling, emphasizing that it should be understood as based 
upon elementary particle/hole excitations of the Mott state. Sec. 0] describes the 
dual vortex description of this conventional integral filling Mott transition, including 
a pedagogical description of hamiltonian boson-vortex duality. Sec. [5] applies the 
dual formulation to non-integral Mott transitions, and describes the origin of the 
main results summarized in this introduction. Finally, Sec. details some specific 
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predictions of the theory for / = 1/2, which provides an example and view to 
"deconfined criticality" , and for / = 1/3, which does not. 



2.1. Physics and model 

The simplest systems that can exhibit a Mott transition between states with 
extended and localized carriers are those composed of interacting lattice bosons, the 
extended states being superfluids/superconductors, and the Mott states being gen- 
erally charge ordered away from integer boson fillings. An exciting experimental 
development of late is the direct realization of such models of cold atomic bosons 
confined to an optical lattice. These will likely turn out to provide the cleanest 
experimental tests of theoretical approaches to this simplest Mott problem. Concep- 
tually, a symmetry-equivalent problem arises if fermions (e.g. electrons) are strongly 
bound into bosonic Cooper pairs. This is not clearly realized in any simple electronic 
material. One may expect, however, that a bosonic theory of this sort will properly 
describe the universal low energy properties of an electronic system in the vicinity 
of a superconductor to Mott insulator transition, provided that both the supercon- 
ductor and insulator are "bosonic". More precisely, we believe it clearly applies 
when the superconductor should be clean (strictly speaking, superclean) and with a 
full gap (e.g. s-wave) to unpaired quasiparticle excitations, and the Mott insulator 
should exhibit only charge (e.g. not spin) order. Similar theories likely apply more 
generally to other superconductor to Mott insulator transitions, but are beyond the 
scope of this paper. 

We concentrate now on purely bosonic Mott transitions. We use a "rotor" rep- 
resentation for the bosons, with phase operators and conjugate number operators 
hi where i runs over the sites of the direct lattice. These operators obey the com- 
mutation relation 



One may consider a variety of geometries, but we will focus here upon the simple 
square lattice. 

A simple boson hamiltonian in the class of interest has the structure 



where t represents a nearest-neighbor boson hopping amplitude, and U an on-site 
boson repulsion. We denote by A a the discrete lattice gradient in the a direction. 
Additional more complicated off-site terms are included in the unspecified contribu- 
tion to the hamiltonian, TC' (but see below). For now we require only that it respect 
the symmetries of the underlying lattice, boson conservation, and locality (e.g. the 
amplitude for hopping between far away sites should decay sufficiently rapidly with 
distance) . 

Eq. (|2-2[) may be studied in the grand canonical ensemble, i.e. at fixed chemical 
potential (2Un). However, one may equally well consider instead the canonical 
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ensemble with fixed boson number (average filling /). We will typically do the 
latter, except in Sees. 12.212.31 and the first part of Sec. |3J where we work at fixed 
chemical potential. 

2.2. Mott states at integral filling 

Neglecting terms in TC', the zero 
temperature phase diagram of TC is well- 
known. 10 ) It takes the schematic form 
in Fig. ^ For t/U <C 1, the system 
is in a Mott insulating ground state, 
with (hi) = N, the integer nearest to 
n, on every site. This phase persists 
inside the "lobes" drawn in the figure. 
There is a gap to the lowest-lying ex- 
cited states, which may be thought of as 
single extra/missing bosons (which delo- 
calize into plane- waves) . For large t/U, Fig 
the ground state is a superfluid (SF in 
the figure), with (e*) = <F sf ^ 0, and 
the density / = (hi) varies smoothly 
with parameters in an unquantized fash- 
ion. There is no excitation gap, and the 
lowest-lying excitations are acoustic "phonons" 
of the broken U(l) symmetry of the superfluid. 




t/u 



1. Schematic phase diagram of boson ro- 
tor model, with on-site interactions only. 
The shaded regions indicate where a there 
is a large near-degeneracy of states with 
different boson densities, and the system is 
highly susceptible to off-site interactions. 

or "phasons", the Goldstone modes 



2.3. Mott states at non-integral filling 

We now return to the shaded regions of the phase diagram in Fig.^ where states 
with different boson density are nearly degenerate. Indeed, in the simple model with 
TC' = 0, for n = iV+1/2, states with any average density between iV and N+l are de- 
generate. For t/U = 0, the eigenvalue of hi = N or hi = N + 1 can be independently 
chosen on each site. The omitted terms in TC' will then clearly determine the nature 
of the ground states appearing in the shaded region. Generally, Mott insulating 
states appear at rational fractional fillings, / = p/q, with p, q relatively prime. For 
q > 1, these are boson "crystals" or charge density waves. Mott states with increas- 
ing q are expected to require longer-range interactions in TC' for their stabilization. 



For example, in the vicinity of n = 
N + 1/2, we can adopt a pseudo-spin 
description, with Sf = hi — N — 1/2 = 
±1/2. In the limit of U -> oo (or t/U <. 
1), one can then replace 

n ^^tJ2(s?S7 + SrSf), (2-3) 

(ij) 

with (ij) indicating the sum is taken 




Fig. 2. Schematic phase diagram with off-site 
interactions. Some representative Mott in- 
sulating states with (fii) — TV + 1/2 are 
shown. 
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over nearest-neighbor sites, and the conventional definitions of spin-1/2 raising and 
lowering operators. This term tends to split the degeneracy of boson states in favor 
of a delocalized superfluid. Terms in TC' can compete with this term for small t, 
splitting the degeneracy instead in favor of a non-trivial Mott state. For instance, 
one may take 

H> = J z SIS* + £ SIS] (S+S-S+Sf + h.c.) + • • • (24) 

Here ((ij)) indicates a sum over next-nearest-neighbors, and the □ a sum over four 
site plaquettes (with sites i,j,k,l labelled clockwise around the plaquette). For 
t <C min{ J z , J' z , K}, a variety of Mott states occur in such a model, including "solids" 
and "valence bond solids" as shown in Fig. |2J The figure also shows a schematic phase 
diagram in which Mott states with q = 2 are stable. More complex phases (with 
q > 2) have been observed in similar models in the literature. 23 -* 

§3. "Integral" Mott-Superfiuid transitions: LGW theory as Bose 

condensation 

It is well-known that the quantum phase transitions across the phase boundaries 
in Fig. ^ are described by LGW theories. 10 ) On symmetry grounds, this is under- 
standable, since the assumptions underlying Landau's theory are valid. In particular, 
the integer filling Mott states, being unique, and having the full symmetry of the 
underlying hamiltonian and only short-range correlations, may be regarded as truly 
"disordered". The superfluid is the "ordered" phase, its symmetry group (the lat- 
tice space group) being a subgroup of the symmetry group of the Mott state (the 
direct product of the lattice space group and the U(l) boson number conservation 
symmetry). A LGW-style expansion of the effective action in terms of the superfluid 
order parameter \P and its derivatives indeed seems to describe the quantum critical 
points. 

It is instructive, however, to understand these transitions more physically. We 
will demonstrate that the physical content of the LGW theory is Bose-Einstein con- 
densation of "particle" and/or "hole" excitations of the Mott state. To that end, let 
us think about the elementary excitations in that phase. For U/t ^> 1, the ground 
state is simply 

\GS)=H(b\) N \0), (3-1) 

i 

where bi = e~ l '^ i is the boson annihilation operator in the rotor formulation, and 
|0) is the state of "no bosons" (zero rotor angular momentum), nj|0) =0 (note that 
eigenstates of hi exist with negative as well as non-negative integer eigenvalues in 
the rotor formulation). In the same limit, the lowest excited states are "particles" 
and "holes" with one extra or missing boson, 

\Pi) = b\\GS), (3-2) 
\hi) = bi\GS). (3-3) 
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For t/U = strictly, the particle and hole energies (relative to the ground state) are 

4°) = 2U(N + \~n), £f = 2U(n + ±-N). (3-4) 

For < t/U <C 1, these states will develop dispersion. By considering the first order 
splitting of the degenerate manifold of particle or hole states (degeneracy associated 
with the site of the particle or hole) , one obtains 

E p/h = - 2t(cos k x + cos k y ) ps A p/h + |— , (3-5) 

where we have Taylor expanded around the minimum at k = 0, giving m = l/{2t). 
The minimum energy excitation is then a particle or hole, for n > N or n < N, 
respectively. The excitation gap, A p / h , is 

A p/h = 2U(~T(n-N))-4t, (3-6) 

to this order. A simple-minded extrapolation of this formula to larger t suggests 
that the excitation gap for particles or holes will vanish once t/U > {\ — \n — 
N\)/2. Though this 0(1) value is well beyond the region of validity of the expansion 
(except for \fi — N\ ps ^, where it is invalid for other reasons, to be discussed below), 
it does suggest a simple physical picture of the transition to the superfluid, as a 
"condensation" of these particles or holes. For n > N, the particles condense, while 
for n < N the holes do. Precisely at n = N, both condense simultaneously. This 
latter case automatically applies in the canonical ensemble if the filling / is fixed to 
be integral. 

We will now show that such particle/hole condensation is the physical content 
of the LGW theory. While it is possible to derive a field theory of this condensation 
from TC, we instead just write it down based on this simple physical picture. We 
model the particle and hole excitations by fields p(x, r), /i(x, r) respectively, in the 
imaginary time (r) path integral. The weight in the path integral is, as usual, the 
Euclidean action, 

S = [ dirfx y(d T +E p -^-V 2 )p+h\d T +E h -^-V 2 )h-X^h^+ph)+- ■■}. (3-7) 
J 2m 2m 

Here we have included a term A which creates and annihilates particles and holes 
together in pairs, which is expected since this conserves boson number. Microscop- 
ically this term arises from the action of the hopping t on the naive ground state, 
which creates particle-hole pairs on neighboring sites, so A ~ 0(t) (the spatial de- 
pendence is unimportant for the states near k = 0). We have neglected - for brevity 
of presentation - to write a number of higher order terms involving four or more bo- 
son fields, representing interactions between particles and/or holes, and other boson 
number-conserving two-body and higher-body collisional processes. Note that the 
dependence upon t/U in Eq. (|3-7|) arises primarily through implicit dependence of 
E p /h- 
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Eq. (|3-7|) can be systematically analyzed by diagonalizing the quadratic form 
involving p, h. The qualitative features can be understood even more simply by 
considering simple limits. First, suppose n > N, so E p < Eh- Upon increasing 
t/U, then, at some point E p — > while Eh > 0. Then the holes are still gapped, 
and the h field can be integrated out order by order in A and higher order terms, 
encountering no vanishing energy denominators. One has then an effective action 
for the remaining particle field p, 



S= did 2 x \p\d T + E P - ^-V 2 )p + u(p^pf + ■■■' 
J 1 2m 



(3-8) 



with E p ~ E p — 0(X 2 /Eh) a renormalized particle gap. This is exactly the well- 
known effective action for the superfluid-Mott transition, which has an LGW form 
with the "order parameter" being the particle field p. The single d T derivative is 
well-understood and present because particle annihilation becomes particle creation 
under time-reversal, and leads to so-called z = 2 behavior of the QCP. 

The case of n < N can be understood similarly with the roles of particles and 
holes interchanged. The remaining case of n = N (the tips of the lobes in Fig. P) 
is slightly different. In that case, both particles and holes become gapless together 
as t/U is increased. Thus we may expect a "relativistic" theory in which particles 
and holes appear in a sense as particles and antiparticles. This is seen by changing 
variables to the linear combinations 

!P=-L(p + /»t) S = -^(p-/»t). (3-9) 

Writing E p = Eh = Eq, one finds that (taking A > without loss of generality) 
the quadratic form for & becomes unstable before that of S, and one can therefore 
integrate out S to obtain 

S = f dufx \&{-dl - —V 2 + rW + u(^V) 2 + ■■■], (3-10) 
J 2m J 

with r ~ Eq — A. Apart from a trivial space-time anisotropy, Eq. ()3-l(J|) is exactly the 
classical LGW free energy of a three-dimensional normal-super fluid transition. Thus 
the two LGW actions, Eqs. (|3-8I3T()|) . usually denoted as particle-hole asymmetric 
and symmetric theories, respectively, are indeed physically equivalent to particle or 
particle+hole condensation. 



§4. "Integral" Mott-Superfiuid transitions: vortex condensation theory 

As discussed in the introduction, a quantum field theory description of the QCP 
would appear to require only some pointlike "particles" described by the quantum 
fields (as creation/annihilation operators). On the Mott side of the transition, parti- 
cle and/or hole excitations provide these, and their condensation precisely coincides 
with LGW theory. On the superfluid side, there is another quite different pointlike 
"particle" excitation: a vortex or anti-vortex. This is specific to two-dimensional 
superfluids, since a vortex becomes a line defect in three dimensions. 
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It is natural then to try to describe the Superfluid-Mott transition, coming from 
the superfmid side by a field theory for vortices. Provided time-reversal symmetry 
T is unbroken, one should expect such a field theory to be relativistic, since vortices 
and anti- vortices are interchanged by T (note that this is independent of the presence 
or absence of particle/hole symmetry in the boson system). One may worry that 
a vortex is a non-local object, with a power-law tail of superflow surrounding it 
extending to infinity. Perhaps this invalidates its use as a particle? 

4.1. Duality 

This worry is resolved by duality, which is a rigorous mathematical mapping 
of the original rotor boson model to one of vortices. We will see that all the non- 
locality of the vortex is taken into account by a non-compact U(l) gauge Geld. 
The dual formulation is mathematically analogous to a lattice U{\) Higgs theory of 
particle physics, or a lattice classical three-dimensional Ginzburg-Landau theory for 
a superconductor (charged superfluid). 

Boson-vortex duality is extensively covered in the literature. It is however often 
a confusing topic for students and senior researchers alike. This is because certain 
stages of the duality mapping are usually performed in a rather non-rigorous fash- 
ion, involving so-called "Villain potentials", fudging with the time-continuum limit 
in the path integral, and other similar manipulations. In reality this sloppiness is 
unimportant, since the usefulness of duality does not lie in the quantitative analysis 
of specific microscopic models, but in its application to universal phenomena. Nev- 
ertheless, it seems anathema to those trained in the more rigorous solid state physics 
tradition rather than the effective field theory approach originating from statistical 
mechanics. 

For this reason we will present here what to our mind is the simplest and most 
rigorous possible version of U(l) duality, performed at the hamiltonian level. This 
proceeds in several steps. The first, completely rigorous and explicit step, contains 
the essence of the duality mapping. It is simply a change of variables from boson 
number and phase fii, <pi living on the direct lattice sites, to new "electric field" and 
"vector potential" variables, E aa ,A aa , living on the links of the dual lattice (We use 
a, b, c to label sites of the dual lattice, and a, f3 = 1, 2 to label links of the dual square 
lattice in the x,y directions, respectively): 

E aa = ^Mh, (4-1) 

Geometrically, the oriented electric field on a link of the dual lattice is taken to equal 
l/(27r) times the phase difference between the phase <pi to immediately to the right 
(using the orientation of the field) of the dual link and the phase <f>j immediately 
to the left of this dual link. The electric field E aa is thus a periodic variable, with 
E aa *-> E aa + 1- The dual vector potential A aa is a discrete field, i.e. has eigenvalues 
of 2ir times integers. It is implicitly defined by Eq. (|4-2f) so that its lattice curl - 
the counter-clockwise circulation around a dual plaquette - is 27T times the boson 
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number inside this dual plaquette. The A aa variables have a "gauge" redundancy: 
the vector potential can be shifted by any 2ttx integer gradient, A aa — > A aa + A a x a 
(with Xa £ 2ttZ) without changing hi. To achieve a one-to-one mapping of A aa to 
hi, one can "fix a gauge" in numerous ways, e.g. for Coulomb gauge, A a A aa = 0. 

It is straightforward to check (by comparing the commutator of the left-hand 
side of Eq. 1)4- 1|) and the right-hand side of Eq. (|4-2|) with the commutator of the 
right-hand side of Eq. (|4T|) and the left-hand side of Eq. (|4-2j) ). that the dual electric 
and vector potential variables are canonically-conjugate: 

[Aaa, E h p] = i5 ab 5 a/3 . (4-3) 

Furthermore, the expression in Eq. (|4T|) implies a constraint: 

A a E aa e Z, (4-4) 

i.e. the lattice divergence of the dual electric field is an integer (a non-zero value 
being "allowed" due to the periodicity of the (pi and E aa variables). Physically, this 
integer can be identified with the vorticity. This can be seen by considering the line 
sum of A a (pi around some area. In an abuse of the continuum notation, 

j> • dl = 2tt j> E ■ dh = 2tt J d 2 x A a E aa , (4-5) 

so that the phase winding around some area, in units of 2tt, just counts the total 
divergence of E aa inside this area, i.e. the net dual "charge" inside this area. 
With full rigor,*) the hamiltonian can be rewritten as 

H = -t cos 2 -*E aa + 77^2 Yl {{^A a Ap)i - 2irn) 2 . (4-6) 

aa i 

Already without further manipulation, Eq. I|4-6|) appears very similar to an elec- 
tromagnetic hamiltonian (oc E 2 + B 2 ). It is, however, inconvenient to work with 
since A aa is a discrete, 2ir x integer- valued field. To understand how to remedy this 
deficiency, it is instructive to express the first term in Eq. (|4-6|) in the A aa basis. 
Since the electric field is conjugate to the vector potential, this cosine is just a "shift" 
operator for the A aa field. Hence we can write 

-tcos27rE aa = -t^2 ([Aaa + 2vr) (A aa \ + \A aa )(A aa + 2vr|). (4-7) 

One can think of the basis for the Hilbert space of a single link of the dual lattice 
as consisting of 2ir x integer-spaced "sites" along a "line" j4 aQ! -space. The original 
boson-hopping term on a link of the direct lattice is just a tight-binding hopping 
hamiltonian for the "line" associated to the dual link crossing this direct link. 

The analogy to a tight-binding model thereby suggests a means of removing the 
discrete constraint on A aa . We can do this by simply replacing the tight-binding 

*' Some small additional care needs to be taken for finite systems with periodic boundary con- 
ditions, in which case additional c- number terms should be added to the left-hand side of Eq. 114-211 . 
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model by a continuum one of a particle in a periodic potential, strategically chosen 
so that the lowest band is (arbitrarily) well-approximated by a tight-binding band, 
well-separated from higher energy states. This amounts to replacing 

-t cos 2irE aa -> —El a - i cos A aa , (4-8) 

and allowing A aa to take arbitrary continuous real values. The tight-binding limit 
is recovered for tn 3> 1. One can then adjust t, k to achieve the desired tight- 
binding matrix element t between nearly-localized levels in neighboring wells of the t 
potential. We do not do this explicitly here, since we will not use it in the following. 
It can, however, be achieved, by taking i, k _1 3> U, to arbitary desired accuracy. 
We are left with the hamiltonian 

H = Y1 ^E 2 aa -t cos A aa +L>^((e Q/3 Z\ Q ^ /3 ) i -2vrn) 2 , (4-9) 

aa L - 1 i 

with U = . This must be supplemented by the commutation relations, Eq. (|4-3j) . 
and the constraint, Eq. (|4-4|) . 

Eq. (|4-9j) is clearly a !7(1) gauge theory, but the presence of "vortex" variables 
is not immediately apparent. In fact, they are implicit in the constraint, Eq. (|4-4|) . 
which can be regarded as Gauss' law (in the dual electromagnetic analogy). Sites 
with non-zero A a E aa thus correspond to dual charges - physical vorticies. Though 
no explicit vortex variables appear, they are un-necessary: the locations of all (dual) 
charges can be determined from the electric field lines. 

It is nevertheless useful (and conventional) to introduce redundant vortex vari- 
ables. In particular, we introduce an auxiliary Hilbert space of integer rotor states, 
eigenvectors of N a G Z, and conjugate variables 6 a , with [N a ,0b] = i5 a b as usual. 
These states are introduced only to redundantly label the vortex positions. So we 
require 

A a E aa = N a . (4-10) 

For this to be consistent with the hamiltonian action in the expanded Hilbert space, 
we must ensure that the — tcosA a ^ term, which changes the electric divergence on 
neighboring sites, also increment the new rotor variables. This is accomplished by 
modifying A aa — > A aa — A a a . Making this shift, one sees that the hamiltonian in 
Eq. (|4-9|) is completely equivalent to the more conventional form 



1 ^2 



2k 



E afl -tcos(A a 6 a - A afl ) 



+ ^^((e^^^) 4 -2vrn) 2 , (4-11) 



combined with the constraint in Eq. (|4-10|) . 

4.2. Phases in the dual formulation, for integer f 

For integer n, we may remove the dual "background flux" 2nn in the last term 
in Eq. 1)4-11(1 by the shift A a 2 — > A a 2 + 2irfa x (with (a x ,a y ) the coordinate of site 
a). Here we are guaranteed / = n by particle- hole symmetry of the on-site rotor 
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model. More generally, working in the canonical ensemble with fixed filling /, this 
shift by construction takes account exactly of the average dual flux in A aa - the 
fluctuations around this shifted around are zero by this choice. One thus has just a 
theory of lattice electromagnetism coupled to a "charged scalar" field tp a = e~ lda . It 
may equally well be thought of as a Ginzburg-Landau-like theory of a dual lattice 
"superconductor" with "pair field" ip a - 

From these analogies, one expects two phases. In the gauge theory language, 
there is a "Coulomb phase" or "dielectric", in which the charged particle is gapped 
and can be integrated out. A simple limit of this phase is obtained by taking t = 0. 
In this limit, N a is a constant of motion on each dual lattice site. Because of the 
Gauss' law constraint, Eq. (|4-1U|) . the ground state is clearly obtained for N a = 0, i.e. 
no vortices present. Individual vortices can be introduced anywhere in the system 
and in this limit have no dynamics, but due to the constraint, an energy cost which 
is logarithmic in system size (due to the dual electric field lines decaying as 1/r far 
from the vortex). Clearly this is the physical superfluid phase in the direct language. 
As in any superfluid, we expect a phason or Goldstone mode. This corresponds to 
the gapless linearly dispersing transverse photon of the dual electromagnetism (there 
is only a single polarization in two spatial dimensions). Going back to the duality 
mapping, it is straightforward to see that small t corresponds to a large "tight- 
binding bandwidth" for A aa , i.e. a large direct boson hopping amplitude, where we 
indeed expect a superfluid state. 

The other phase is a "Higgs phase" or dual "superconductor", in which the 
V> a particle is condensed. In this "Higgs" phase the photon is gapped, and indeed 
there is an energy gap to all excitations. Moreover, there is no broken symmetry: 
the condensate amplitude (tp a ) is not gauge invariant and the vortex condensation 
does not itself represent an order parameter for any broken symmetry. This then 
corresponds to the featureless Mott insulating state in the direct picture. 

4.3. Continuum field theory (for integer f) 

For integer /, a rather naive continuum limit is possible, due to the ability to 
transform away the background "flux" in Eq. I|4-llj) . It is instructive, paralleling the 
logic used above to derive the LGW theory coming from the Mott state, however, to 
think more physically in terms of the vortex excitations. We imagine coming from 
the superfluid state, by increasing t starting from a small value. At first pass, we 
will neglect the fluctuations of A aa , E aa . 

For t = 0, the ground state |0) clearly has N a = 0, and corresponds to the vortex 
"vacuum". The lowest excited single vortex states correspond to N a = ±1 on an 
(arbitrary) single site a, which we denote 



These states are elementary in that any vortex number configuration can be built 
from them by superposition. Each such state has a logarithmically-divergent energy, 
since the Gauss' law constraint, Eq. (|4-1U[) . requires E aa ~ 1/r far from the vortex: 
this is just the usual logarithmic vortex energy coming from the long-range superflow. 




(4-12) 
(4-13) 




: 0, 
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We will return to this point below. 

First order degenerate perturbation theory in t gives an effective single-particle 
tight-binding model amongst the |a±) states, and hence some plane- wave eigenstates 
with energy dispersion. One may now follow the same logic as in Sec. |3J with the 
vortex/anti- vortex states playing the analogous roles to the particle/hole excitations 
in the Mott phase. Time-reversal symmetry guarantees that these have the same 
excitation gap. Recognizing that these states carry the dual gauge charge, and 
allowing for fluctuations in the A aa ,E aa fields, one thereby arrives (following the 
same methodology as in Sec. 121 using the coherent state path integral, and the 
standard path integral Trotterization of electromagnetism) at the continuum action 

Sduai = f drcfx [^(e^xd u A x f + - iA^rf + r \<pf + u\^\ A + •••]. (4-14) 

Here for compactness, we have neglected to write unimportant space-time anisotropies 
(constant scale factors) that appear between space and time derivatives, and spatial 
(electric) and temporal (magnetic) fluxes, writing simply three-dimensional coor- 
dinates (t,x,u) and indices fj,,f, A = 0,1,2. The anisotropies do not lead to any 
additional physics, and are expected to be irrelevant or redundant at the QCP. New 
coupling constants r, u, e 2 have been defined in the continuum action. 

Eq. (|4-14|) is identical to the classical free energy of a three-dimensional Ginzburg- 
Landau theory for a superconductor. The duality between this form and Eq. (|3-1U|) . 
the classical three-dimensional XY free energy, was established in Ref. 12). Because 
the LGW form, Eq. (|3-lUj) . does not involve any gauge field, it is of more practical 
(analytical and numerical) use in this case. Nevertheless, we emphasize that both 
the LGW and dual actions are descriptions of the same critical point. Qualitatively 
the same physics can be extracted from either form, though they will yield different 
quantitative results from approximate treatments. 

It is instructive to review how the measures of charge localization are encoded 
in Eq. (|4-14|) . The boson creation operator, b\, is a 2ir flux insertion operator in the 
dual formulation, i.e. it creates a space-time monopole with d^B^ = 2ir5^ (x, r), 
where the dual space-time magnetic flux = e^\d u A\ is physically the boson 3- 
current. By looking at monopole correlators, therefore, one can discern the presence 
or absence off-diagonal long-range order, or the presence or absence of a gap. The 
superfluid density p s and compressibility k measure the response to external physical 
vector and scalar potentials, respectively. These physical potentials couple to the 
boson 3-current, and hence p s and k are obtained from correlation functions of the 
A^ gauge field. The lesson to be learned here is that all Mott/Superfluid properties 
are encoded in the properties of the gauge field A^, not directly in the dual "order 
parameter" ip (which anyway is not gauge-invariant - See Sec, l4,3.Tl for further related 
discussion). 

A studious reader may wish to try the following illustrative exercise: calculate 
the critical behavior of the superfluid density using the dual form in the "random 
phase approximation". This can be done by introducing an external gauge field, 
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A^ xt , which couples to the physical electromagnetic current, 

S -» S - i Jdrd 2 x A e *e^xd u A x . (4-15) 

Neglecting the \ip\ 4 term, integrate out 99 at quadratic level in A^, then integrate 
out A^ itself, to obtain the coefficient of (PrA^) 2 (Pt is a transverse projection 
operator), which is proportional to p s . One may compare the manner in which this 
vanishes with f with the mean-field LGW prediction for how p s vanishes with r in 
Eq. (ETTOl) . 

4.3.1. Logarithmic vortex potential and bound states 

A careful reader may be troubled by the notion of using a vortex as an elementary 
excitation, since, in the superfluid state, it has a logarithmically divergent energy 
(in the system size). A partial answer to this question is that any neutral collection 
of equal numbers of vortices and antivortices has finite energy. However, it is clear 
that a vortex and anti-vortex attract each other, and will form a bound state, with 
infinite (logarithmically) binding energy. On increasing t within the superfluid state, 
one may suspect that such a vortex/anti- vortex pair excitation will condense before 
individual vortices do. 

In fact, more careful thought is needed. In the "relativistic" theory of Eq. (|4-14|) 
(or the lattice hamiltonian, Eq. (|4-11|0 vortices and anti- vortices are not separately 
conserved. Indeed, the action of t on the ground state creates them in pairs on 
nearest-neighbor links (compare to the A term in Eq. (|3-7|0 . Thus the number of 
such "bound states" is not a conserved quantity, and consequently there is no sharp 
phase transition associated to their "condensation" - any "creation operator" for 
these bound states has a non-vanishing expectation value for all t > 0, unless it is 
finely tuned. Thinking spectrally, were such a single neutral bound state to approach 
zero energy upon increasing t, it would be expected to have an avoided crossing with 
the ground state, because of the non-zero matrix element between the naive ground 
state and the state with one excited bound pair. 

More generally, Eqs. (|4-14l4-llj) have no internal symmetries, only space-time 
symmetries and "gauge invariance" which is not a symmetry at all but embodies 
a dynamical constraint. Phase transitions in these models are then expected to be 
either associated with the development of a Higgs mass for the gauge field, or by 
breaking of spatial lattice symmetries (the latter not being expected to occur for 

fez). 

Nevertheless, it is a reasonable physical question to ask why it is that a de- 
scription in terms of elementary individual vortex excitations is appropriate when 
these objects are always infinitely strongly bound on the superfluid side of the transi- 
tion. The answer is that, as the Superfluid-Mott QCP is approached, the interaction 
between vortices is becoming progressively more and more "screened" by virtual fluc- 
tuations of vortex/anti- vortex pairs. This can be seen, for instance, from the fact 
that the superfluid density vanishes at the QCP. Alternatively, just from scaling, at 
the critical point, the interaction between two static external dual "test charges" 
separated a distance r behaves like 1/r, not logarithmically. A screened interaction 
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prevails on scales smaller than the correlation length, which diverges as the QCP is 
approached. 

§5. Non-integral Mott-Superfluid transitions: vortex theory 

Having understood the direct and dual formulations of the integral Superfluid- 
Mott transition, we turn to the non-integral case. We will concentrate on rational 
mean fillings, / = p/q, with p, q relatively prime, and q > 1 but not too large. 
We assume that sufficient interactions are present to stabilize a "crystalline" Mott 
insulating state at such a density. We must then analyze Eq. (|4T1|) . with n ~ p/q 

m) 

If we attempt to approach the problem from an LGW perspective, there is a 
fundamental difficulty: neither the Mott insulator or Superfluid are "disordered" 
phases, the former breaking space-group symmetries of the lattice, and the latter 
breaking the U(l) boson conservation symmetry. Landau theory can only describe 
the transition from one of these states to an even more symmetry-broken phase, 
which then seems to require an intermediate phase between the two. 

Taking a more physical point of view, we can search for pointlike excitations in 
either phase that could provide the degrees of freedom for a critical quantum field 
theory. Unfortunately, on the Mott side, the nature of the elementary excitations 
would seem to be very specific to the particular Mott state under consideration. 
All the Mott phases presumably include extra/missing boson "particles", as well as 
domain wall and other discrete topological excitations particular to the precise type 
of boson density order in the ground state. The extra/missing boson "particles" could 
provide the variables for some QCP, but one expects condensation of such particles 
not to disrupt the density-wave order of the Mott state, leading to a "supersolid" 
rather than a true superfluid. As such states appear to be even more exotic than 
superfluids and Mott insulators, we will not explore that possibility further. A theory 
based on topological excitations of particular Mott states is possible, but much more 
limited in scope (see Ref. 1), Sec. ?? for a discussion). 

The most general approach seems to be from the superfluid. This is quite ap- 
pealing insofar as, unlikely the panoply of possible Mott states, the superfluid is 
unique, and has a symmetry which is independent of the boson filling. Thus the 
dual approach, based upon vortex degrees of freedom, generalizes directly to arbi- 
trary filling factors. Moreover, we will see that the variety of distinct Mott states 
arises naturally from this description. 



*' Strictly speaking, in the superfluid, the density is not exactly equal to n, except when n is an 
integer or half integer in the microscopic rotor model. What is important is not n but the density 
/ in the superfluid adjacent to the Mott state (within which / is independent of n for a range of n, 
it being incompressible). The actual value of n should be adjusted to achieve / = p/q. We tacitly 
assume this is done below. 
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5.1. Vortex projective symmetry group 

As for integral filling, we shift the dual vector potential to include the average 
dual flux, A aa — ► A aa + A aQ , with 

A a2 = 2vr/a x , A al = 0, (5-1) 

corresponding to Landau gauge. Note that the specific gauge choice (and indeed 
any gauge choice) breaks naive spatial lattice symmetries. Nevertheless, because 
e a pA a A a p represents a uniform gauge flux, this is an artifact of gauge fixing. Clearly, 
since A aa is a dynamical variable, it is possible simply to undo the effect of any space 
group operation on A aa by an appropriate shift of A aa . Because the gauge-invariant 
flux is not changed by such an operation, this shift of A aa must be pure gauge, 
i.e. a gradient A a x a , for some scalar \a- As a consequence, we can compose this 
shift with a pure gauge transformation - of both the vortex field a (or equivalently 
-0a) and the gauge field - by the phase factor Xa in order to undo the shift of A aa . 
By this reasoning, it is clear that for every operation in the space-group, there is 
a corresponding transformation consisting of a naive space-group transformation 
and a gauge rotation of the vortex fields without any transformation of the gauge 
field. This transformation is almost unique, up to a global (i.e. a-independent) U(l) 
phase rotation of the vortex fields. The global phase arbitrariness is allowed because 
a uniform gauge rotation has zero gradient and hence does not shift A acx . 

In this way, by making some arbitrary phase choice for each group element, 
one associates a transformation with each operation in the space group. These new 
transformations generate some new group, which obeys the original group multipli- 
cation table up to phase factors, which in general cannot be removed. This is called 
a projective representation of the space group, and in a slight abuse of notation, we 
also call this new group a Projective Symmetry Group, or PSG. The proper subgroup 
(i.e. not including reflections) of the PSG on the square lattice is generated by the 
three operations associated with x and y translations, and a ir/2 rotation (which we 
choose about a dual lattice site), which can be chosen as: 

T y : if){a x , a y ) — > i[)(a x ,a y - 1) 



Tjdual 
K tt/2 



tp(a x , a y ) -> ip{a x - 1, a y )uj ay 

t/j(a x ,a y ) -> ip(a y ,-a x )uj a * a y, (5-2) 



with to = e 2m f . Specific algebraic relations follow from these definitions, notably 

T x T y = u)T y T Xl (5'3) 

which informs us that, unlike in the original space group, the operations associated 
with translations in the PSG do not commute. The translational subgroup of the 
PSG is well-studied, and Eq. (|5-3|) is known as the magnetic translation algebra. The 
physical meaning of Eq. ()5-3(l is, as explained in the introduction, that the vortex 
acquires a dual Aharonov-Bohm phase of 2irf upon encircling a site of the direct 
lattice, containing on average / bosons. 
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The remaining algebraic relations between generators are unchanged from the 
original space group: 



5.2. Vortex multiplets 

We again consider now the nature of the elementary vortex excitations, starting 
from Eq. ()4-llj) . As in Sec. 14.31 at leading order in t, we must solve the tight-binding 
model arising from degenerate perturbation theory amongst the |a±) states, this time 
in the presence of the mean gauge field, Eq. (|5-lj) . This is the Hofstadter problem, 
which is well-known generally to have a extremely complex "butterfly" spectrum. 24 ) 
We are, fortunately, interested only in the lowest-energy states at the bottom of 
the lowest Hofstadter band. As for any particles, these must appear in multiplets 
comprising an irreducible representation of the symmetry group of the hamiltonian. 
In this case, this group is the PSG. 

It is straightforward to see that, for a filling/dual flux with denominator q > 1, 
all representations of the PSG are at least (/-dimensional. This follows because one 
may choose to, say, diagonalize T y . However, by Eq. ()5-3j) . acting with T x on a state 
multiplies its eigenvalue of T y by u>, hence this must be a linearly independent state 
to the initial one. This can be done repeatedly q — 1 times, until the q th time, one 
arrives back at a state with the initial eigenvalue of T y . Thus the PSG connects 
states of q different values of the quasimomentum. Since we used only Eq. (j5-3j) . 
this conclusion is clearly true for the translational subgroup of the PSG alone, so 
representations of the full PSG can only be equal or larger. 

It turns out, from explicit solution of the tight-binding model, that a (/-dimensional 
irreproducible representation (irrep) does exist for the full PSG, and moreover, this 
is the lowest energy vortex multiplet for the case of interest. Details of this construc- 
tion can be found in Ref. 1). The result is that q vortex fields, ipi, £ = . . . q — 1 may 
be defined, as relativistic field operators (i.e. analogously to Eq. ()3-9j) . as superpo- 
sitions of "particle" (vortex) creation and "anti-particle" (anti-vortex) annihilation 
operators) for each member of the multiplet. Under translations, 



Here, and henceforth, the arithmetic of all indices of the ipi fields is carried out 
modulo (/, e.g. (p q = tp$. Under the rotation, 




(5-4) 



T x : ipe — > (fe+i 



(5-5) 




(5-6) 



which is a Fourier transform in the space of the q fields. The full (improper) PSG is 
generated by including also reflections (defined here about x and y axes of the dual 
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lattice) , 

If al : n - Vi 

If al : ^ - ^ . (5-7) 

A continuum field theory can be now constructed from the g-dimensional vortex 
multiplet. The effective action should be invariant under all the physical symmetries, 
with vortex fields transforming under the PSG and the dual U(l) gauge symmetry. 
Restoring fluctuations of the gauge field, the critical action constructed in this man- 
ner, generalizing Eq. (|4-14|) to non-integer filling, has the form given in Eq. IJ1-2JI in 
the introduction. 

5.3. Order parameters 

It is important to consider the observables in the problem. For integer filling, 
we discussed that ip is not itself a true "order parameter" because it is not gauge 
invariant. Physical quantities (e.g. superfluid density, compressibility, off-diagonal 
long-range-order) in that case are related just to properties of the dual gauge field. 
This is because local, gauge-invariant combinations of the single vortex field such as 
|</?| 2 etc. (related to the vortex/anti- vortex bound states discussed in Sec. I4.3.l"]l are 
scalars under all symmetries of the hamiltonian. 

The situation is dramatically different for q > 1. While the above physical 
quantities are still related to properties of the dual gauge field, in general, gauge- 
invariant bilinears of the form y'lw are n °t scalars under the spatial symmetries. 
They can thus serve as order parameters for various types of symmetry breaking. 
Group theoretically, the direct product of an irrep of the PSG and its conjugate can 
be decomposed into a sum of true irreps of the ordinary space group (not projective 
representations). Happily, the basic components of these irreps can be constructed 
in generality. 1 - 1 With the definition of Eq. IJ1-4JI of the introduction, straightforward 
manipulations show that p(Q) transforms exactly as expected for a Fourier compo- 
nent of a scalar "density" with wavevector Q. Specifically, it is easy to verify that 
Pmn = P-m.-n, and from Eqs. (|5-5|) and (|5-6[) that the space group operations act 
on p mn just as expected for a density wave order parameter 

Ti _ — m 

x • Pmn ' Pmn 

Ty . Pmn * & Pmn 

-^7r/2 : Pmn > P—n,m- (5'8) 

§6. Examples 

While the set of different p mn can describe a variety of different density wave 
orders in the Mott state, the number of such ordering patterns is limited, so the 
vortex theory actually constrains the nature of Mott insulating density ordering 
occuring in the neighboring of a continuous transition to a superfluid. These orders 
can be determined from a mean-field analysis of the effective action S, Eq. (|4-14|) . 
We give examples for q = 2,3. Further examples can be found in Ref. 1). 
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6.1. Half-filling and deconfined criticality 

The case of q = 2, corresponding to bosons at half-filling, is particularly inter- 
esting. One has Ant = £4 + 0(ip e ), with 

r (\ |2 1 I 1 2\ 2 , 701 / * * \2 ,\ 

£ 4 = — (Ivol +lvi| J + - • l 6 ' 1 ) 

It is convenient to make the change of variables 

C0 + C1 



■ Co - Ci / c oA 

<pi = -t ■=—. (6-2) 



V2 
)_ 

V2 



The action in Eq. ()6-lj) reduces to 



A = ^(|Co| 2 + |Ci| 2 ) 2 -^(lCo| 2 -|Ci| 2 ) 2 . (6-3) 

The result in Eq. (|6-3|) is identical to that found in earlier studies 15 )' 16 ) of the q = 2 
case. 

Minimizing the action implied by Eq. ()6-3j) . it is evident that for 701 < there 
is a one parameter family of gauge-invariant solutions in which the relative phase of 
Co and Ci remains undetermined. As shown in earlier work, 15 ) this phase is pinned 
at specific values only by an 8 th order term proportional to ~ (CoC* ) 4 + c - c - 

The mean-field analysis finds three phases: First, phase (A): an ordinary charge 
density wave (CDW) at wavevector (w, tt). The other two states are VBS states in 
which all the sites of the direct lattice remain equivalent, and the VBS order appears 
in the (B) columnar dimer or (C) plaquette pattern. The phases (A), (B), and (C) 
appear in the upper-left, upper-right, and lower-right corner of the inset in Fig. |5J 
Note that the "site-centered stripe" phase in the lower-left of the inset, though it 
certainly can occur in lattice boson models, does not occur according to the vortex 
theory in the vicinity of the transition to the superfluid. The saddle point values of 
the fields associated with these states are: 



(A) 
(B) 
(C) 



Co + , Ci = or Co = , Ci / 0. 
Co = e m " /2 Ci + 0. 

Co = e i(«+l/2W2^ ^ Q) (64) 



where n is any integer. 
6.1.1. Deconfined criticality 

This particular example, for 701 < 0, has recently been understood in much 
more detail. It turns out that this case, describing the transition from a superfluid 
to a columnar or plaquette VBS state, can be understood from a complementary 
point of view as a theory of fractional "half" -boson excitations. It is too involved to 
fully explore this in detail in this paper, but we will at least uncover these fractional 
excitations in the dual theory. 
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For 701 < 0, the mean-field solutions in the Mott state have both Q = |C|e 4 ^, 
with |£| constant. Supposing slowly varying the effective phase-only action is 



S. 



eff 



A ^ 2 + 2#( e i"* dt ' A > 



(6-5) 



with p' s = 2|£| 2 . Consider a fixed 'vortex' (independent of r) in one - say i?o ~~ °f the 
2 phase fields, centered at the origin r = 0. One has the spatial gradient V^o = 4>/r, 
while i?i = (here </> = (—y,x)/r is the tangential unit vector) . Clearly, the action 
is minimized for tangential A = Acp. Far from the 'vortex' core, the Maxwell term 
{(-iivxdvA\) 2 is negligible, so one need minimize only the first term in Eq. 16-51 The 
corresponding Lagrange density at a distance r from the origin is thus 



C'v = % [(1/r - Af + A 2 



(6-6) 



Minimizing this over A, one finds A = l/(2r). Integrating this to find the flux gives 



fx 



dr 



7T 



(6-7) 



Since the physical charge is just this dual flux divided by 2ir, the 'vortex' in (q indeed 
carries fractional boson charge 1/2. 

Note that, because of the small 8 th order term locking the phases i?o and i?i 
together, such a 'vortex' in just one of these fields costs a divergent energy (actually 
linear in system size). This indicates these fractional particles are "confined" in the 
Mott state, like quarks in quantum chromo dynamics. This confinement, however, 
becomes weaker and weaker as the SF-Mott QCP is approached, due to the smallness 
of the 8 th order term in this limit. In this sense - and in others discussed in Refs. 19 ) 
- this is a "deconfined QCP". 

A naive extension of this argument would suggest the presence of charge 1 jq for 
general q. It turns out that this generalization is not so simple, and while it can be 
made in some cases, there are considerable restrictions involved - see Ref. 1). 

6.2. q = 3 

Now there are 3 ipi fields, and the quartic potential in Eq. H6-1JI is replaced by 



+ 



7oo 
4 

7oi 



2\2 



Oo| +|</>i| +M 2 ) 
(<Po<P*<f% + PiV&o + <P2<Po ( P-<PoVi + c - c - 



- 2|< y 2 | 2 |^i| 2 - 2|< / ji| 2 |<£ 2 | 2 - 2|(£2| 2 |<A)| 2 ) • 

The results of a mean-field analysis for this potential are shown in Fig 01 
The states have stripe order, one 

along the diagonals, and the other along 

the principle axes of the square lattice. " 1 " 1 



(6.8) 
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Both states are 6-fold degenerate, and the characteristic saddle point values of the 
fields are 

(A) : tpo / , (px = ip 2 = or 

e i4«V3^ 2 = e i2nir/ 3lpi = ^ ^ q 

(B) : ^ = W = e ±2 " /3 ^2 

and permutations. (6'9) 
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